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Abstract - For piecewise expanding one-dimensional maps without periodic turning points 
we prove that isolated eigenvalues of small (random) perturbations of these maps are close to 
isolated eigenvalues of the unperturbed system. (Here "eigenvalue" means eigenvalue of the 
q^ | corresponding Perron-Frobenius operator acting on the space of functions of bounded vari- 

ation.) This result applies e.g. to the approximation of the system by a finite state Markov 
chain and generalizes Ulam's conjecture about the approximation of the SBR invariant mea- 
sure of such a map. We provide several simple examples showing that for maps with periodic 
turning points and for general multidimensional smooth hyperbolic maps isolated eigenvalues 
are typically unstable under random perturbations. Our main tool in the ID case is a spe- 



ed 

Jl ■ 

cial technique for "interchanging" the map and the perturbation, developed in our previous 



paper |6|, combined with a compactness argument. 



1 Introduction 

We discuss stochastic stability in the following general framework: A discrete time dynamical 
system is a pair (/, X), where X C H is a bounded phase space (say X = [0, l] d ) and 
/ : X — ► X is a nonsingular map, iterations of which define trajectories of the system. 
Nonsingular means that m(f~ 1 A) > for any measurable set iCI with positive Lebesgue 
measure m(A) > 0. 

Consider now small random perturbations of the discrete time dynamical system. Roughly 
speaking this means that, when we apply / to a point x £ X, rather than choosing the exact 
value of fx we choose in a random way, in accordance with some distribution, a point from 
the ball B £ (fx) (i.e. with centre at the point fx and radius e). 

Definition 1.1 Let Q e (x,A) be a family of transition probabilities and f : X — » X a map. 
We denote by f e the Markov process on the phase space X defined by the transition probabil- 
ities Q £ (fx,A) and call f £ a random perturbation of f. 

*On leave from Russian Academy of Sciences, Inst, for Information Transmission Problems, B. Karetnij 
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Our main assumptions on the perturbations are the following: 

Q e (x,A)=0 if dist(», A) > e (locality), (1.1) 

\\Q £ h-h\\ 1 < d(Q £ ) -var(/i) -► as e ^ (smallness), (1.2) 

vav(Q E h) < vax(h) + C \\h\\ x (regularity) (1.3) 

for any function h of bounded variation (var(/i) < oo). The parameter e here plays the role of a 
"magnitude" of the perturbation. These assumptions are satisfied for a broad class of random 
perturbations, among them convolutions with absolutely continuous transition probabilities, 
bistochastic absolutely continuous perturbations, singular perturbations of point mass type, 
deterministic perturbations by chaotic maps close to identity, and Ulam type perturbations 
(see, for example, @ for details). 

Stochastic stability of the Sinai-Bowen-Ruelle (SBR) measure of a dynamical system 
discussed under the same assumptions in Q may be considered as a weak kind of stability, 
because other statistical characteristics may be unstable. In this paper we further explore 
the problem of stochastic stability and study the stability of the spectrum of the Perron- 
Frobenius operator Pj (PF-spectrum for brevity), considered as an operator on the Banach 
space (BV, ||-|| ) of functions of bounded variation. This operator describes the dynamics 
of densities under the action of the map /. Our main stability result provides rather general 
sufficient conditions for stochastic stability of isolated eigenvalues of the spectrum in the one- 
dimensional piecewise expanding case, and counterexamples below show that these conditions 
can hardly be relaxed. Moreover we demonstrate that arbitrarily small random perturbations 
and especially Ulam type perturbations (see below) of a generic multidimensional hyperbolic 
map (having a stable manifold) can completely change the PF-spectrum of the system. 

One of the main motivations for the present work was to study the stability of the PF- 
spectrum in Ulam's construction of a finite state Markov chain approximation [1^, |j of 
chaotic dynamics. The idea of the construction is to take a finite partition {Aj} of the phase 
space with bounded volume ratios and to approximate the action of the map / by the Markov 
chain with transition probabilities 

lA.nr^-l 

Pij := — fAi — • 

This construction can be considered as a special type of small random perturbations, where 
the transition operator satisfies all above assumptions. The convergence of invariant measures 
of these finite Markov chains to the SBR measure of the approximating dynamical system 
under this construction was proved in [fbf l for piecewise expanding (PE) maps with large 
enough expanding constants, and for general PE maps in || ||] (see also a review there). 
For numerical applications of this method see e.g. ||] . Since Ulam's construction provides a 
rather general approach for numerical modeling of chaotic dynamics, the question of stability 
of the PF-spectrum and the study of nonconvolution type random perturbations becomes 
important not only from a purely theoretical point of view but also from a practical one. 

The operator Pj- preserves integrals, thus its leading eigenvalue is equal to 1 and the most 
important feature of its spectrum S(Py) is the modulus of the second largest spectral value 
r2 := sup{|t| : r G £(Pj), |r| < 1}, which characterizes the rate of convergence to the SBR 
measure. The corresponding value for the perturbed operator is denoted by r^. 



The basic idea here is the following. Recall that \\h\\ = var(/i) + 11^11!- In || we proved 
that the transition operator P £ = Q £ Pf of the randomly perturbed map satisfies the uniform 
Lasota-Yorke type inequality 



P?h 



<a-\\h\\ +C-||/ l || 1 (1.4) 



BV 



for h £ BV, some fixed integer N, a E (0, 1) and C > independent of e. This yields at once 
the existence of /V-invariant densities h F with \\h F \\ < 4= L — , such that (11.21) forces each weak 



limit /i* = lim e ^o h e to be an invariant density for /. (In fact, as | \h £ \ | is uniformly bounded, 
/i* is not only a weak limit, but ||/i* — h £ \\ l — > as e — ► 0.) If the expanding constant A 
is larger than 2, ( |1.4| ) was proved with N = 1 in various settings, see e.g. [J10], |3L Q and 
references therein. For quite a while it was supposed that the extension of this inequality to 
cases with A G (1, 2] is only a technical problem. However, the counterexample constructed in 
[Q] shows that the situation is not so simple. After this counterexample it became clear, that 
the main problem is the possible existence of periodic turning points, i.e. points where the 
derivative of the map / is not well defined. Namely, under the action of random perturbations 
"traps" or "absorbing sets" can appear near these periodic turning points, which leads to the 
appearance of new localized ergodic components in the perturbed system. 

Definition 1.2 Let X = [0, 1]. A map f : X — > X is piecewise C 2 if there exists a partition 
of X into disjoint intervals {Xj}, such that /Iciosfx ) * s a C 2 -diffeomorphism (of the closed 
interval Clos(Xj) to its image). Its expanding constant is defined as 

A/:= .inf \f'(x)\. 

j, x£Xj 

A piecewise C 2 map is called piecewise expanding (PE), if Xrk > 1 for some iterate f k . 

Definition 1.3 The image of a measure /i under the action of a map f is the measure 
ffi defined by f ^{A) = n(f~ 1 A) for any measurable set A. By f e \i we mean the measure 
f e fi(A) = f Q £ (fx,A)dfj 1 (x). A measure n is f (f £ ) -invariant, if ffj, = fi (f £ fi = fi). /i is 
called smooth, if it has a density with respect to Lebesgue measure. 

Definition 1.4 A turning point of a map f is a point, where the derivative of the map is 
not well defined. We denote the set of turning points by TP, and the set of periodic turning 
points by PTP. 

Recall that the essential spectral radius of the operator Pj is the smallest nonnegative 
number 8 for which elements of the spectrum S(Pj) outside of the disk of radius 9 centered 
at the origin are isolated eigenvalues of finite multiplicity. In [11] it was shown that for a PE 
map 6 = linin^oo (7 A7„ . Therefore / is PE if and only if 8 < 1. From now on we fix some 
numbers 8' , 8 arbitrarily close to 8 with 8 < 8' < 8 < 1. 

Theorem 1.1 Let f be a piecewise expanding map with PTP = 0, and let perturbations Q £ 
satisfy (\1. 1\ ), !\l.t\ ) and (1.5). Let r be an accumulation point of eigenvalues r £ > 8 > 8 of 



the perturbed operators P £ := Q £ Pf for e — > 0. Then there are a sequence k —>■ oo and a 
function h £ BV such that r £k — > r and \\h £k — h\\ 1 — > as k — ► oo ; where P//i = rh and 



f£ k h £k — r £k h £k . 



In other words, any accumulation point of the eigenvalues of the perturbed operators 
lying outside of the disk containing the essential spectrum of the original operator is an iso- 
lated eigenvalue of Pt . This generalizes the corresponding part of the results for convolution 
type perturbations of smooth expanding maps in |2| . In the case of Ulam type approxima- 
tions based on Markov partitions of smooth hyperbolic maps our theorem is complemented 
by the following result |J: Each isolated eigenvalue is a limit point of eigenvalues of the 
corresponding Ulam operators. 

It is worth to remark that the rate of convergence of the eigenfunctions for the eigenvalue 
1 is 0(e|loge|) whereas we have no rates of convergence for other spectral quantities. In 
particular, in some cases numerical experiments show very slow convergence of the r £ to true 
eigenvalues. 

Technically the proof of this theorem is based on the following proposition, being a con- 
sequence of several technical results obtained in || . 

Proposition 1.1 Let f be a piecewise expanding map with PTP = 0, and let perturbations 
Q £ satisfy ( \1. 1\ ), \1.^\ ) and ( \1. $) . Then there exist constants C, Eq such that for some finite 
N 

var(Pf h) < e N var(h) + C | 



for any e 6 (0,£o) an d any function h £ BV. 

The paper is organized as follows. In Section 2 we recall some necessary definitions and 
prove our main stability result. Section 3 is devoted to the analysis of various situations 
when the spectrum is not stable with respect to random perturbations. Especially important 
among these situations are the spectrum collapse in the absence of isolated eigenvalues and 
the instability of the spectrum of multidimensional hyperbolic maps due to the presence of 
the stable foliation of the map. Finally we discuss a possible generalization of the notion of 
the spectrum by means of zero- noise limit of the spectra of randomly perturbed systems. 

2 Proof of the stability result 

Recall that the variation of a function over h : X = [0, 1] — ► Hi is defined as var(/i) := 
sup{ f x c/)'hdx}, where the supremum is taken over all functions <fi £ C 1 (1R) with compact 
support, H^Hoo < 1 and H^'H^ < oo. Notice that as X = [0,1] is bounded, \\h\\ 1 < ^var(/i) 
for all h £ BV, and if I C X is an interval, then var(/i • 1/) < var(/i). Indeed, our setting 
means that var(/i) is the variation of h ■ lx over H. 

Given N and (3 > as defined above we refine the partition Z into intervals of mono- 
tonicity of the map / by adding further points to TP in such a way that 



var 



(X^ z f lz ), var((A / | Z / | ' z )- 1 )</?, (2.1) 



not introducing new TP that are mapped to other TP. 
Define P 1} P 2 : BV -> BV by 

P 1 h = Q £ P f (h-l x \ Y ), P 2 h = Q £ P f (h-l Y ) 

where Y is a neighbourhood of TP scaling linearly with respect to e (see 



Proposition 2.1 /^, Proposition 3.1] Suppose that there are constants Ci,C 2 > and a G 
(0, 1) such that 

var(Pf /t) < da k vai{h) + C 2 \\h\\ x for all k G 7L + (2.2) 

and that there is some N G 2Z+ such that 

P 2 P\P 2 = forallk = l,...,N. (2.3) 

Then 

var((Q e P)"/i) < ( iV(iV 2 +1) C 1 3 + d) / • var(/ l ) 

+(l + C 1 + C 1 2 )C 2 ^±^-||/ i || 1 . 

The assumptions fl2.2| ) and ([D|) were verified in Q for PE maps with PTP = and for 
sufficiently small e without paying attention to a particularly sharp estimate of the constant 
a. For the purposes of our present paper we need to show that a can be chosen as a = 8', 
i.e. close to the essential spectral radius 6. 

Since PTP = 0, there is some no such that P^ ° = so that inequality (|2.2| ) can be 
satisfied for any positive a with suitable constants C\,C 2 depending only on / and Q. For 
Pi the relevant estimate is given in || Proposition 3.8]. The value a = (f) 1 ' given in the 
statement of the proposition is not as sharp as the corresponding proof permits. In fact, the 
effictive estimate derived in the proof is 

var(P 1 iV /i) 

< (1 + ^) 2N (X^ + l -Np\- f N )^{h) + {Xf N + l -Np\f)C N ■ \\h\l, 

for each h G BV. [] 

Proof of Proposition |1.1| . It remains to choose the constant (3 > so small that 

(1 + IpfN {Q N + l _ N p X -N ) < {df) N < ~ 6 N 



to obtain the statement of Proposition 1.1. 



Proof of Theorem |l.l| . Suppose that P £ h £ = r £ h £ for some h £ G BV with 1 1 /i^ 1 1 ^ = 1 and 
|^e| > 0' . Then there is a constant S > depending only on 9 and on the constants in the 
Lasota-Yorke type inequality such that var(/i e ) < S. Indeed, by Proposition |Ll| 

var(P £ JV /i £ ) <^var(/i £ ) + C||/i e || 1 . 

Since P^h £ = r^h £ , we obtain the following estimate: 

i f e\ N C 

var(/i e ) = var(p 7 P e Ar /i e ) < ( — j var(/i e ) + -^ \\h £ \\ 1 , 



1 Observe that the brackets in the corresponding formula in flq] are set slightly wrong, which had no further 
effect on the subsequent proofs in that paper. 



or 

C 

var(/i e ) < — =— := S < oo, 

v £ > ~ r N _ QN 

since r £ > 9. In particular also var(P^/i e ) < S and var(Pj/i) < S, where h = ho. 

Let r be any accumulation point of the r £ for e — » 0. Then there are a sequence e^ — » 
and a function /i 6 51/ such that r £fc — ► r and \\h £k — h\\ l — > as k — ► oo. In particular 
var(/i) < S. Using the assumption ( |1.2| ) it follows that 

\\p f h-p £k h £k \\ 1 <\\(p f -p Ek )h\\ l + wp^h-h^w, 

<d(Q £fc )-S 1 +||/ i -/ l£ J| 1 ^0 



as e — > and 



Hence P//t = rh. 



\rh — r £ ,h £ .\\, < \r\ ■ \\h — h £l IL + |r — r e , I • ||/i e . Ik — > . 
i c k c kiii — ii ii ^kiii i & fc i ii ^kiii 



3 Instability of the PF-spectrum for generic maps 

3.1 Instability of the essential spectrum 

The above estimates ensure stochastic stability of the isolated eigenvalues of Pj in the sense 
that spectral points outside the essential spectrum cannot arise "from nothing" under per- 
turbations. We have no examples of maps / without PTP and perturbations Q £ where an 
isolated eigenvalue of Pj is not approximated by eigenvalues of the P £ . However, if the 
only isolated eigenvalue is 1, the essential spectrum may collapse, which leads to spectrum 
localization (?2 = 0). The simplest example when this phenomenon takes place is the case 
of Ulam perturbations of the well known dyadic map x — ► 2x (mod 1) , when the number of 
intervals in the Ulam partition is equal to 2 n . If this number is not an integer power of 2 
there will be a stable eigenvalue fi = \j1 (on the boundary of the essential spectrum) with 
the eigenvector (x — 1/2). The following statement generalizes this observation. 

Theorem 3.1 Let f : X — » X C 1R be piecewise linear with respect to a partition Z and 
assume that f(I) = X for each I G Z. Denote by Zn the partition of X into domains of 
linearity of f . Then the spectrum of the Ulam operator Pz N constructed with respect to the 
partition Z^ consists of only and 1. 

Sketch of proof: Let /ibea piecewise constant function on Zn such that / h(x) dx = 0. 
Then P% h = 0, which proves that r2 = 0. 

3.2 Instability in the periodic turning points case 

In this more general case we cannot proceed as above in spite of the fact that the first 
eigenvector is stable under some additional restrictions (see jg] for details). The reason is that 
the proof of this statement in || uses slightly weaker estimates of the rate of convergence, 
which do not ensure our Proposition [Ll]. In this section we want to argue that isolated 

6 



eigenvalues of the spectrum (except 1) might be unstable indeed. To show this we modify 
the W-map (described in detail in ||) such that locally around the fixed turning point c 
it remains the same, while no point "outside" is mapped into a small neighborhood of the 
point c. The modified W-map is shown in the Figure [y. Observe that this map is Markov. 
Consider the following random perturbation: x — > x + e with probability 1 — 5 and x — ► x with 
probability 5. We apply this random perturbation only in the neighborhood of the point c 
consisting of two neighbouring intervals of monotonicity. For 5 = there is a e-small ergodic 
component I £ around the point c. For small enough e, 5 > the escape rate out of the set 
I £ is of order 5. By the construction of the map, the probability to hit into the interval I e 
starting outside of it is equal to zero. Therefore the transition operator corresponding to the 
randomly perturbed system must have an eigenvalue of order 1 — 5. Thus the spectral gap 
of the perturbed operator vanishes as 5 — > 0. 

Example 1. Let us discuss now in more detail Ulam perturbations of the modified W- 
map. Let x be a local coordinate such that the fixed turning point be 0, and denote the 
Ulam interval containing this point by / = [—a, b]. Define A = \f! |. Then the probability to 
remain in the inteval I = [—a, b] in Ulam's approximation is equal to 

a/A + min{a, Xb} IX a + min{a, Xb} 

R a > °) = ; = w t~\ ■ 

v ' a + b X(a + b) 

Let a = Aft. Then 

p(AM) = x(xbTV) = XTT = 1 - ATT" 

On the other hand, if a < b then the image of any other Ulam interval does not intersect 
with the interval I. This shows that 2/ (A + 1) < 1 belongs to the spectrum of the perturbed 
operator, whereas the special form of the modified W-map implies that the spectrum of P/, 
except for the eigenvalue 1, is contained in the circle with radius 1/A < 2/(A + 1). Observe 
that in this case the perturbation cannot make the spectral gap arbitrarily small, but r<i 
might be unstable. In any case this proves the appearance of an eigenvalue of the perturbed 
operator in the "spectral gap" of the unperturbed one. 

One might argue that still it is possible that under the action of well behaving random 
perturbations of convolution type the influence of just a few periodic turning points may 
not matter. To show that this is not the case, consider the following symmetric random 
perturbation: 

{x — e with probability q 
x with probability 1 — 2q 

x + e with probability q, 

with < q <C 1. We apply this pertutbation to the same modified W-map for A = 2 (notice, 
that in this case r2 = 1/2), and show that the corresponding transition oprator Q £ P has an 
eigenvalue 1 — 2g > 1/A. Since no point outside of a small neighborhood of the fixed point 
can hit into this neighborhood, it is enough (as in our previous examples) to study the escape 
rate from this neighborhood. If A = 2 then locally the behaviour of the randomly perturbed 
system is completely described by the following random walk model on 7L: 

if x> 
otherwise, 





Figure 1: The modified W-map. 



-1 with probability q 

with probability 1 

1 with probability q. 



2q 



Clearly, if \x\ > 1 then the trajectory of this point will never return to zero. Consider the 
part of the transition matrix corresponding to the points —1,0 and 1. It can be written as 




Thus this matrix has an eigenvalue 1 — 2q, which proves our statement. 



3.3 Instability for a generic multidimensional hyperbolic map 

Stability of spectral properties becomes a much more delicate problem in the multidimensional 
case. Traditionally, to define this spectrum one considers the Perron-Frobenius operator 
for the expanding map defined on unstable manifolds induced by the original map (see, for 
, [?]] ) . Another way to calculate the isolated eigenvalues is to study the so called 



example, [16 



weighted dynamical (-function of a map, which counts periodic points of the map weighted by 
Jacobians in the unstable direction. Zeros and poles of the ("-function correspond to isolated 
eigenvalues of the map (see, for example, |l]] and references therein). 

Our spectral stability results can be generalized for finite systems of weakly coupled ID PE 
maps (see || for definitions). On the other hand, even for a more general multidimensional 
PE map our construction does not work, since there is no good control over coefficients of a 
Lasota-Yorke type inequality in this case, contrary to the ID case. 

Example 2. Consider a smooth hyperbolic map / : 1R — > H . Let be a hyperbolic 
fixed point of the map, and let the horizontal direction be locally unstable with the expanding 
constant \ u > 1, while the vertical direction be contracting with A S <1. We consider Ulam 
partitions into equal squares rotated by the angle ir/2 with respect to the coordinate axes. 



/I \\ *u 


V 



Figure 2: One element of the rotated Ulam partition and its preimage in the 2D hyperbolic 
case. 



One element of the partition together with its preimage is shown in Figure |2| Straightforward 
calculations show that the probability to remain in the considered Ulam square is of order 
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One can calculate this probability exactly also for the case of finite values of \ s , for example 
p = 2/3 for A s = 1/2. This is not sufficient to prove that fi > T2, but it indicates that the 
limit behaviour of the approximation might differ from that of the original map. A slightly 



less striking example of this type was discussed in [13], where it was claimed that the worst 
situation is when the angle between the contracting and expanding directions is small. 

The following numerical example shows that the untypical instability of the essential 
spectrum due to the presence of periodic turning points in the one-dimensional case becomes 
typical for multidimensional maps. Near a periodic point of a multidimensional hyperbolic 
map stable and unstable foliations are coming arbitrarily close one to another. Therefore 
an arbitrarily small (random) perturbation can mix them (similarly to the situation near 
periodic turning points (see 0])). We study an example as simple as possible to demonstrate 
that this type of behavior is generic. 

Example 3. Consider the well known "cat" map, which is the simplest example of a 
smooth two dimensional hyperbolic map. This is a map from the unit torus X = [0, 1] x [0, 1] 
into itself defined by (x, y) \— ► (x + y (mod 1) , x + 2y (mod 1) ). We consider two partitions 
of X into equal squares. First we simply divide horizontal and vertical axes into n equal 
intevals, whose products give a partition into N = n 2 squares, which we call the standard 
partition. There is a one to one correspondence of these squares and pairs of integers (i = 
nx,j = ny), where (x,y) is the pair of coordinates of the lower left corner of a square. Here 
i,j £ {0, 1, . . . ,n — 1}. Simple calculation gives the following transition probabilities for the 
corresponding Markov chain whose elements are numbered as jn + i + 1 (see also Figure |3|): 




Figure 3: Image of an element of the "standard" partition by the "cat" map. 





i + j,i + 2j 


i + j, i + 2j + 1 


i + j + l,i + 2j + l 


i+j + l,i + 2j + 2 


i,3 


1/4 


1/4 


1/4 


1/4 



Moduli of the "second" eigenvalues (r2) of the transition matrices (n 2 x n 2 ), and their 
multiplicities (in parenthesis) are shown in the following table: 



n 


r-i 


n 


T2 


n 


T2 


n 


T2 


n 


T2 


2 

7 

12 

17 


0.0000(3 ) 
0.4886(16) 
0.3783(48) 
0.4335(24) 


3 

8 
13 
18 


0.3536(8 ) 
0.4454(24) 
0.4835(28) 
0.5957(24) 


4 

9 

14 

19 


0.0000(15) 
0.3847(24) 
0.4886(16) 
0.5387(36) 


5 
10 
15 
20 


0.3299(20) 
0.4275(12) 
0.4045(16) 
0.4275(12) 


6 
11 
16 

21 


0.3536(8 ) 
0.5161(20) 
0.4454(24) 
0.5357(32) 



Compare this with the inverse to the largest eigenvalue of our linear map 1/A = 2/(3 + 
\/5) w 0.38204, which (see, for example, J?]]) is the correct value of the "second" eigenvalue 
in this case. 

To show that even this is not the worst case we consider also another partition, namely 
the standard partition shifted by l/(2n) (in both directions). Similarly to the previous case, 
we associate the square centered at (x,y) with the pair of integers (i = nx,j = ny). Here 
i,j £ {0,1,..., n — 1}. The transition probabilities for the corresponding Markov chain 
(whose elements are numbered as jn + i + 1) are shown in the following table: 





i+j,i + 2j 


i + j, i + 2j + 1 


i + j + l,i + 2j + l 


i + j,i + 2j-l 


i + j-l,i + 2j-l 


i,j 


1/2 


1/8 


1/8 


1/8 


1/8 



Moduli of the "second" eigenvalues (r2) of the transition matrices (n 2 x n 2 ), and their 
multiplicities (in parenthesis) are 



n 


T2 


n 


T2 


n 


T2 


n 


T2 


n 


T2 


2 


0.3968(3 ) 


3 


0.3953(8) 


4 


0.4543(12) 


5 


0.4029(20) 


6 


0.4443(24) 


7 


0.5577(16) 


8 


0.4940(24) 


9 


0.5038(24) 


10 


0.4754(12) 


11 


0.6203(20) 


12 


0.4567(48) 


13 


0.5371(28) 


14 


0.5577(16) 


15 


0.5495(16) 


16 


0.4940(24) 


17 


0.5864(36) 


18 


0.6733(24) 


19 


0.5976(36) 


20 


0.4902(24) 


21 


0.5838(32) 
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In this case all "second" eigenvalues are greater in modulus than 1/A. Looking at these 
two tables the structure of limit points of the eigenvalues seems not quite clear for both of 
the considered families of the partitions, and one might argue that for large enough n the 
corresponding r2 may converge to 1/A. However, we have numerical evidence that for the 
standard partition for n = Ik there is an eigenvalue 0.4886, while for the shifted standard 
partition for n = 8k there is an eigenvalue 0.4940. The following general statement justifies 
this prediction and provides us with the precise description of the structure of the spectrum 
for the case of a linear automorphism preserving integer points. 

Theorem 3.2 Let f : H d — ► JR d be a linear map such that f{7L d ) = 7L d , and let the map f := 
f (mod 1) be defined on the d-dimensional unit torus. Denote by P n the matrix corresponding 
to Ulam's approximation of the map f constructed according to the partition of the unit torus 
into n d equal cubes. Then r G S(P^ n ) for any positive integer k whenever r E E(P n ). 

The proof of this result is based on the fact that due to the selfsimilar structure of Ulam's 
approximation in this case both the matrix P^ n and the eigenvector e^n corresponding to the 
eigenvalue r consist of repeated blocks of the matrix P n and the eigenvector e n respectively. 

Therefore if for some n we obtain numerically a "bad" value for the "second" eigenvalue, 
it will still be present for large enough multiples of n. 

In recent papers || it was proposed to use Ulam's procedure based on a finite Markov 
partition to estimate r<i. This claim was justified in these papers for ID smooth expanding 
maps and 2D Anosov automorphisms. In practice, the usefulness of this approach is limited 
by the observation that usually such partitions can be found only numerically, and as we 
shall show a small error here may lead to even worse accuracy of the eigenvalues compared 
to a non Markov partition. 

Now we are in a position to answer the question why the spectral gap in the above 
numerics differs significantly from theoretical predictions. To have a simple model for the 
analytical study to start with, consider a family of 2D maps from the unit square into itself, 
defined as follows: / 7 (x, y) := (2x (mod 1) , 7(2/ — c) + c (mod 1) ), 7 > 0, < c < 1. For each 
value of 7 the map / 7 is a direct product of two ID maps. Thus, for 7 > 1 one can prove 
that the PF-spectrum of this map is just the set of all pairwise products of elements of the 
spectra of the involved ID maps. For 7 < 1, however, only the first map contributes to the 
spectrum. Observe that for 7 < 1 the invariant measure is concentrated on the attracting 
fiber r := {(x,c) : < x < 1}, and the spectrum is the PF-spectrum of the piecewise 
expanding map on T. The dependence of r2 on the parameter 7 is shown by thick lines in 
Figure ||. Observe the discontinuous behaviour when the parameter 7 crosses the value 1. 
By dots we indicate the rate of correlations decay with respect to the Lebesgue measure in 
this system. Observe that these two graphs differ only for < 7 < 1, i.e. when the stable 
foliation is present. 

This simple example demonstrates the main difference between expanding and hyperbolic 
maps, because the "traditional" spectrum in the latter case does not take into account the 
behaviour of the system along the stable foliation. 

3.4 Random perturbations of contractive maps. 

The above examples show that in order to understand how small random perturbations change 
the behaviour of a hyperbolic system one has to study their influence on a pure contractive 
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Figure 4: 7-2(7) f° r the direct product map. 

map. Let / 7 , c (^) := l{x — c) + c, < 7, c < 1, be a family of maps from the unit interval [0, 1] 
into itself. If our random perturbation has the transition probability density q(-, •), being a 
BV function of the first variable, then the corresponding transition operator is well defined 
as an operator from BV into itself and one can compute its spectrum. Denote by [x] the 
closest integer to the point x, and let 5 = \[cn] — cn\ G [0, 1/2] be the distance from the fixed 
point c to the closest end-point of the Ulam interval to which it belongs, multiplied by n. 
The transition matrix P n is lower triangular in this case. Therefore the diagonal entries of 
the matrix are just its eigenvalues. A simple calculation gives the following representation 
for T2 as a function of 5: 



r 2 (S) 




1) 



if S = and 7 > 

if0<<5< T ^ 

otherwise. 



This result shows the following. First, fi sensitively depends on the distance to the closest 
end-point of the Ulam interval it belongs to: r2(0) = 2 — ^ f2(0+) = 1, while f2(l/2) = 
(3 — l/7)/2 (provided 7 > 1/2). Second, when the fixed point lies very close to the boundary 
of one of the Ulam intervals, the estimate is the worst, which yelds a very bad accuracy if 
one uses an approximation to a Markov partition for Ulam's procedure. 

Let us show that shift-invariant random perturbations may cure this pathology. Suppose 
the random perturbation has a shift-invariant transition probability density q(x, y) = q(y—x), 
q G C 2 . We want to advocate that the PF-spectrum of the perturbed operator in zero noise 
limit in this case is well defined, does not depend on the shape of q(-) and is nontrivial. Let 
us prove this for fi- 

Lemma 3.1 Let q £ C 2 , q(x) = if \x\ > e, and e < max {7c, 7(1 — c)}. Then fi = 7. 



Proof. The random map can be rewritten as x n _|_i = j(x r 
sequence of iid random variables with probability density q(- 



X n +l 



l n (xi 



c)+c+ ( 7 n ^i + 7 n ' 2 6 + 



- c) + c + £ n , where (£ n ) is a 
Therefore 



Let £("' := c+J2k=o l k ^n~k- Since the £& are iid, the sequence {£} n ') converges in distribution 
to a random variable ^'°°'. Then the random variables x n converge in distribution to ^°°' 
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as n — ► oo and fi corresponds to the rate of this convergence in the following sense: ^°°' 
can be rewritten as ^( 00 )=^ n (^( 00 ) — c) + ^ n ' , where = means equality in distribution and 
£(°°) is a copy of £(°°) which is independent of the £( n >. Denote by q n ,qoo, h n , <p n and tp n 
the densities of the random variables £, , 6 j x n> 7 n (xi ~~ c ) and 7 n (£'°°' — c) respectively. 
Then /i n +i = q n * 4>n and ^ = q n * ip n , and the supports of n and ip n are of order 7™. 
Observing that for any h £ (7 1 



|/i(x + <5) - /i(x)| cfe = (|(J| + o((5))var(/t), 
we conclude that for large n 

ll<?n-<7oo|li = / \q n {x) - q n (x - y)\ip n (y) dy dx < 0(~f n )v&r(q n ), 

var(,„ - fc ) = / |A(.) - i(, - ,)*,(,)! * d* < OCVWi), 

and similarly 

||<?n - /in+illi < 0(7 n )var(g n ), var(g n - /i n+ i) < 0(7 n )var(g n ). 

As £' n ) = £ n + 7^' n_1 \ there is some probability density h such that q n = q-kh. Therefore 
var(g ra ) < var(g) and var(g^) < var(g') so that 

\\K ~ ^00^ < 0(7™) • (var(g) + varfaO). 

It is easily seen that one can choose the initial density h\ such that this order of convergence 
is attained (e.g. take h\ close to a 5- function) . This yields the claim of the lemma. g 
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